** MAD and Taboo Replication File **

* This file contains replication commands for tables and figures in MAD and Taboo: U.S. Expert Views on Nuclear Weapons"
* This file requires the "estout" Stata package to run. To install type: "ssc install estout" (without quotes) into the Stata command window. 

** Table 1 in the paper reports hypotheses, there is no data to populate the table **

** Table 2 **
* Policymaker Survey Demographics *
* Note: These have been delinked to one another and responses to protect anonymity of respondents *

clear

use "MAD and Taboo Replication Data Final.dta", clear

sum age 
tab sex 
tab minority 
tab military  
tab rank 
sum experience 
tab education  
tab university 
tab thinktank 
tab privateindustry 
tab advocacy 
tab international 
tab political 
tab ideology 

* Replication Command for Note B in Table 2. Note: Denominator used for Note B is the total number answering nuclear familiarity questions (320) *  
* Note: It is necessary to use the separate file that contains the data on university affiliation linked to education *

clear

use "MAD and Taboo University Education Replication.dta", clear

tab education if university == 1

** Replication Command for Footnote 6 **
* Note: It is necessary to use the main data file *

clear

use "MAD and Taboo Replication Data Final.dta", clear

ttest taboo_uk_france, by(taboo_uk_france_early) welch
ci means taboo_uk_france_1
ci means taboo_uk_france_2

ttest taboo_russia_china, by(taboo_russia_china_early) welch
ci means taboo_russia_china_1
ci means taboo_russia_china_2

ttest mad_confidence, by(mad_confidence_early) welch
ci means mad_confidence_1 if mad_confidence_1 != 9
ci means mad_confidence_2 if mad_confidence_2 != 9

ttest deter_nuclear_strikes, by(deter_nuclear_strikes_early) welch
ci means deter_nuclear_strikes_1 if deter_nuclear_strikes_1 != 9
ci means deter_nuclear_strikes_2 if deter_nuclear_strikes_2 != 9

ttest coerce_nuclear_opp, by(coerce_nuclear_opp_early) welch
ci means coerce_nuclear_opp_1 if coerce_nuclear_opp_1 != 9
ci means coerce_nuclear_opp_2 if coerce_nuclear_opp_2 != 9

ttest deter_conv_nuclear_opp, by(deter_conv_nuclear_opp_early) welch
ci means deter_conv_nuclear_opp_1 if deter_conv_nuclear_opp_1 != 9
ci means deter_conv_nuclear_opp_2 if deter_conv_nuclear_opp_2 != 9

ttest coerce_nonnuclear_opp, by(coerce_nonnuclear_opp_early) welch
ci means coerce_nonnuclear_opp_1 if coerce_nonnuclear_opp_1 != 9
ci means coerce_nonnuclear_opp_2 if coerce_nonnuclear_opp_2 != 9

ttest deter_conv_nonnuclear_opp, by(deter_conv_nonnuclear_opp_early) welch
ci means deter_conv_nonnuclear_opp_1 if deter_conv_nonnuclear_opp_1 != 9
ci means deter_conv_nonnuclear_opp_2 if deter_conv_nonnuclear_opp_2 != 9

** Table 3 Demographics **
* Faculty Survey Demographics *

tab sex_ir if deter_nuclear_strikes_ir != . 
tab minority_ir if deter_nuclear_strikes_ir != . 
tab rank_ir if deter_nuclear_strikes_ir != . 
tab ideology_ir if deter_nuclear_strikes_ir != . 

** Figure 1 **
* Confidence in MAD *
* NOTE: Data for 2011 Policymaker survey available at: https://carnrank.nd.edu/files-that-pertain-to-the-above-project/data/" 

tab mad_confidence

** Figure 2 **
* Nuclear Deterrence **

* Policy *

tab deter_nuclear_strikes
tab deter_conv_nuclear_opp
tab deter_conv_nonnuclear_opp

* Faculty *

tab deter_nuclear_strikes_ir
tab deter_conv_nuclear_opp_ir
tab deter_conv_nonnuclear_opp_ir

** Figure 3 ** 
* Nuclear Coercion *

* Policy *

tab coerce_nuclear_opp
tab coerce_nonnuclear_opp

* Faculty *

tab coerce_nuclear_opp_ir
tab coerce_nonnuclear_opp_ir

* Supporting t-test for Figure 3: Policy v. Faculty views on Coercion of Nuclear Opponents *
* NOTE: Joint Coerce Nuclear Opponent drops all 'Don't Knows' *

ttest j_coerce_nuclear_opp, by(policy_survey) welch
ci means coerce_nuclear_opp if coerce_nuclear_opp !=9
ci means coerce_nuclear_opp_ir if coerce_nuclear_opp_ir !=9

** Table 4: Policymaker Ordered Logit Estimates **
* NOTE: To protect respondent anonymity this command cannot be replicated with the files.  The code used is placed here to show the analysis.  *
*eststo M1: ologit deter_nuclear_strikes age sex minority military rank experience education university international political ideology if deter_nuclear_strikes !=9, robust
*eststo M2: ologit deter_conv_nuclear_opp age sex minority military rank experience education university international political ideology if deter_conv_nuclear_opp !=9, robust
*eststo M3: ologit deter_conv_nonnuclear_opp age sex minority military rank experience education university international political ideology if deter_conv_nonnuclear_opp !=9, robust
*eststo M4: ologit coerce_nuclear_opp age sex minority military rank experience education university international political ideology if coerce_nuclear_opp !=9, robust
*eststo M5: ologit coerce_nonnuclear_opp age sex minority military rank experience education university international political ideology if coerce_nonnuclear_opp !=9, robust

*esttab M1 M2 M3 M4 M5, nodepvars nomtitles se(3) pr2 b(3) star(* 0.05 ** 0.01 *** 0.001) title(Table 4: Ordered Logit Estimates of Policymaker Views on Nuclear Deterrence and Coercion) varwidth(50)

** Table 5: IR Faculty Ordered Logit Estimates **

eststo M6: ologit deter_nuclear_strikes_ir sex_ir minority_ir rank_ir ideology_ir if deter_nuclear_strikes_ir !=9, robust
eststo M7: ologit deter_conv_nuclear_opp_ir sex_ir minority_ir rank_ir ideology_ir if deter_conv_nuclear_opp_ir !=9, robust
eststo M8: ologit deter_conv_nonnuclear_opp_ir sex_ir minority_ir rank_ir ideology_ir if deter_conv_nonnuclear_opp_ir !=9, robust
eststo M9: ologit coerce_nuclear_opp_ir sex_ir minority_ir rank_ir ideology_ir if coerce_nuclear_opp_ir !=9, robust
eststo M10: ologit coerce_nonnuclear_opp_ir sex_ir minority_ir rank_ir ideology_ir if coerce_nonnuclear_opp_ir !=9, robust

esttab M6 M7 M8 M9 M10, nodepvars nomtitles se(3) pr2 b(3) star(* 0.05 ** 0.01 *** 0.001) title(Table 5: Ordered Logit Estimates of IR Faculty Views on Nuclear Deterrence and Coercion) varwidth(50)

** Figure 4 **
* Confidence in Taboo *

* Policy *
tab taboo_uk_france
tab taboo_russia_china

* Faculty *

tab taboo_uk_france_ir
tab taboo_russia_china_ir
tab taboo_us_ir

** Figure 5 **
* Difference of Means in Taboo *

* Policy *
ttest taboo_confidence, by(uk_france) welch
ci means taboo_uk_france
ci means taboo_russia_china

* Faculty *
ttest taboo_confidenceufcr_ir, by(uk_france_ir) welch
ci means taboo_uk_france_ir if taboo_uk_france_ir != 9
ci means taboo_russia_china_ir if taboo_russia_china_ir != 9

ttest taboo_confidenceufu_ir, by(us_ir) welch
ci means taboo_us_ir if taboo_us_ir != 9

ttest taboo_confidencecru_ir, by(us_ir) welch

** Footnote 13 Replication **
* Note all variables exclude Don't Know responses *

* Mann-Whitney test *

* Policy *

ranksum taboo_confidence, by(uk_france) 

* Faculty * 

ranksum taboo_confidenceufcr_ir, by(uk_france_ir) 
ranksum taboo_confidenceufu_ir, by(us_ir) 
ranksum taboo_confidencecru_ir, by(us_ir) 

* Nonparametric equality of means test *

* Policy *

median taboo_confidence, by(uk_france) 

* Faculty *

median taboo_confidenceufcr_ir, by(uk_france_ir) 
median taboo_confidenceufu_ir, by(us_ir) 
median taboo_confidencecru_ir, by(us_ir) 

*** Online Appendix Replication Commands ***

* Note: It is necessary to use the correct data files for each figure *

* Figure A.1 Policy - Sex *

clear

use "MAD and Taboo Figure A1 Policy_Sex Replication.dta", clear

tab deter_nuclear_strikes if sex == 1
tab deter_conv_nuclear_opp if sex == 1
tab deter_conv_nonnuclear_opp if sex == 1
tab coerce_nuclear_opp if sex == 1
tab coerce_nonnuclear_opp if sex == 1

tab deter_nuclear_strikes if sex == 0
tab deter_conv_nuclear_opp if sex == 0
tab deter_conv_nonnuclear_opp if sex == 0
tab coerce_nuclear_opp if sex == 0
tab coerce_nonnuclear_opp if sex == 0

* Figure A.2 Faculty - Sex *

clear

use "MAD and Taboo Replication Data Final.dta", clear

tab deter_nuclear_strikes_ir if sex_ir == 1
tab deter_conv_nuclear_opp_ir if sex_ir == 1
tab deter_conv_nonnuclear_opp_ir if sex_ir == 1
tab coerce_nuclear_opp_ir if sex_ir == 1
tab coerce_nonnuclear_opp_ir if sex_ir == 1

tab deter_nuclear_strikes_ir if sex_ir == 0
tab deter_conv_nuclear_opp_ir if sex_ir == 0
tab deter_conv_nonnuclear_opp_ir if sex_ir == 0
tab coerce_nuclear_opp_ir if sex_ir == 0
tab coerce_nonnuclear_opp_ir if sex_ir == 0

* Figure A.3 Policy - Race *
* NOTE: To protect respondent anonymity this command cannot be replicated with the files.  The code used is placed here to show the analysis. *

*tab deter_nuclear_strikes if minority == 1
*tab deter_conv_nuclear_opp if minority == 1
*tab deter_conv_nonnuclear_opp if minority == 1
*tab coerce_nuclear_opp if minority == 1
*tab coerce_nonnuclear_opp if minority == 1

*tab deter_nuclear_strikes if minority == 0
*tab deter_conv_nuclear_opp if minority == 0
*tab deter_conv_nonnuclear_opp if minority == 0
*tab coerce_nuclear_opp if minority == 0
*tab coerce_nonnuclear_opp if minority == 0

* Figure A.4 Faculty - Race *

tab deter_nuclear_strikes_ir if minority_ir == 1
tab deter_conv_nuclear_opp_ir if minority_ir == 1
tab deter_conv_nonnuclear_opp_ir if minority_ir == 1
tab coerce_nuclear_opp_ir if minority_ir == 1
tab coerce_nonnuclear_opp_ir if minority_ir == 1

tab deter_nuclear_strikes_ir if minority_ir == 0
tab deter_conv_nuclear_opp_ir if minority_ir == 0
tab deter_conv_nonnuclear_opp_ir if minority_ir == 0
tab coerce_nuclear_opp_ir if minority_ir == 0
tab coerce_nonnuclear_opp_ir if minority_ir == 0

* Figure A.5 Policy - Ideology *

clear

use "MAD and Taboo Figure A5 Policy_Ideology Replication.dta"

tab deter_nuclear_strikes if ideology3 == 3
tab deter_conv_nuclear_opp if ideology3 == 3
tab deter_conv_nonnuclear_opp if ideology3 == 3
tab coerce_nuclear_opp if ideology3 == 3
tab coerce_nonnuclear_opp if ideology3 == 3

tab deter_nuclear_strikes if ideology3 == 2
tab deter_conv_nuclear_opp if ideology3 == 2
tab deter_conv_nonnuclear_opp if ideology3 == 2
tab coerce_nuclear_opp if ideology3 == 2
tab coerce_nonnuclear_opp if ideology3 == 2

tab deter_nuclear_strikes if ideology3 == 1
tab deter_conv_nuclear_opp if ideology3 == 1
tab deter_conv_nonnuclear_opp if ideology3 == 1
tab coerce_nuclear_opp if ideology3 == 1
tab coerce_nonnuclear_opp if ideology3 == 1

* Figure A.6 Faculty - Ideology *

clear

use "MAD and Taboo Replication Data Final.dta", clear

tab deter_nuclear_strikes_ir if ideology3ir == 3
tab deter_conv_nuclear_opp_ir if ideology3ir == 3
tab deter_conv_nonnuclear_opp_ir if ideology3ir == 3
tab coerce_nuclear_opp_ir if ideology3ir == 3
tab coerce_nonnuclear_opp_ir if ideology3ir == 3

tab deter_nuclear_strikes_ir if ideology3ir == 2
tab deter_conv_nuclear_opp_ir if ideology3ir == 2
tab deter_conv_nonnuclear_opp_ir if ideology3ir == 2
tab coerce_nuclear_opp_ir if ideology3ir == 2
tab coerce_nonnuclear_opp_ir if ideology3ir == 2

tab deter_nuclear_strikes_ir if ideology3ir == 1
tab deter_conv_nuclear_opp_ir if ideology3ir == 1
tab deter_conv_nonnuclear_opp_ir if ideology3ir == 1
tab coerce_nuclear_opp_ir if ideology3ir == 1
tab coerce_nonnuclear_opp_ir if ideology3ir == 1

** Table A.1: Policymaker Familiarity with and Confidence in the Nuclear Taboo **
* Note: There were no Don't Know responses for any category *
* NOTE: To protect respondent anonymity this command cannot be replicated with the files.  The code used is placed here to show the analysis. *
*eststo M11: logit taboo_familiar q1 sex minority military rank experience education university international political ideology if taboo_familiar != ., robust
*eststo M12: ologit taboo_russia_china q1 sex minority military rank experience education university political international ideology, robust
*eststo M13: ologit taboo_uk_france q1 sex minority military rank experience education university political international ideology, robust
*esttab M11 M12 M13, nodepvars nomtitles se(3) pr2 b(3) star(* 0.05 ** 0.01 *** 0.001) title(Table A.1: Logit and Ordered Logit Estimates of Policymaker Views on the Nuclear Taboo) varwidth(50)

* Results reported in Appendix text discussing education and academic affiliation and relationship to familiarity with the nuclear taboo *
* NOTE: To protect respondent anonymity this command cannot be replicated with the files.  THe code used is placed here to show the analysis. *
*tab taboo_familiar if education == 5 & university == 1
*tab taboo_familiar if education != 5 | university != 1

** Table A.3 Results for demographics and taboo views among IR faculty **

eststo M14: ologit taboo_russia_china_ir sex_ir minority_ir rank_ir ideology_ir if deter_nuclear_strikes_ir !=9, robust
eststo M15: ologit taboo_uk_france_ir  sex_ir minority_ir rank_ir ideology_ir if deter_conv_nuclear_opp_ir !=9, robust
eststo M16: ologit taboo_us_ir sex_ir minority_ir rank_ir ideology_ir if deter_conv_nonnuclear_opp_ir !=9, robust
esttab M14 M15 M16, nodepvars nomtitles se(3) pr2 b(3) star(* 0.05 ** 0.01 *** 0.001) title(Table A.2: Ordered Logit Estimates of IR Faculty Views on the Nuclear Taboo) varwidth(50)


